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Using a sample of 1.31 billion J/f) events collected with the BESIII detector at the BEPCII collider, the 
decays J/f) —>■ and J/f) —>■ (/)7r®7r®7r® are investigated. The isospin violating decay J/fi —>■ 
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4>tv° fo{980) with /o(980) — >■ tttt, is observed for the first time. The width of the /o(980) obtained from 
the dipion mass spectrum is found to be much smaller than the world average value. In the 7r'^/o(980) mass 
spectrum, there is evidence of /i (1285) production. By studying the decay J/ip ^ (j>r}', the branching fractions 
of r)' —>■ 7r'*’7r“7r° and T]' —>■ as well as their ratio, are also measured. 

PACS numbers: 13.25.Gv, 14.40.Be 


I. INTRODUCTION 

The nature of the scalar meson /o(980) is a long-standing 
puzzle. It has been interpreted as a qq state, a KK molecule, 
a glueball, and a four-quark state (see the review in Ref. [1]). 
Further insights are expected from studies of /o(980) mix¬ 
ing with the 00(980) [2], evidence for which was found in 
a recent BESIII analysis of J/tjj and Xci decays [3]. BE- 
SIII also observed a large isospin violation in J/ij} radiatively 
decaying into 7r+7r“7r° and involving the interme¬ 

diate decay ?7(1405) —^ 7r°/o(980) [4]. In this study, the 
/o(980) width was found to be 9.5 ± 1.1 MeV/c^. One pro¬ 
posed explanation for this anomalously narrow width and the 
observed large isospin violation, which cannot be caused by 
00(980) — /o(980) mixing, is the triangle singularity mecha¬ 
nism [5, 6]. 

The decays J/'0 —>■ (/)7r+7r“7r° and J/ij) —>■ (()7r°7r°7r° are 
similar to the radiative decays J/ij) —>■ 77r“''7r“7r°/7r°7r°7r° as 
the (j) and 7 share the same spin and parity quantum numbers. 
Any intermediate /o(980) would be noticeable in the tttt mass 
spectra. At the same time, a study of the decay J/ip ^ cpr]' 
would enable a measurement of the branching fractions for 
if —>■ 7r+7r“7r'’ and 7' —>■ 7r°7r°7r°. The recently measured 
B{r]' 37r°) = (3.56 ± 0.40) x lO-^ [4] from a study of the 

decay J /ip ^ ^rj' was found to be nearly Aa higher than the 
previous value (1.73 ± 0.23) x 10“^ from studies of the reac¬ 
tion 7r“p 71(67) [7-9] ^ Additionally, the isospin-violating 
decays 77' —>■ 7r+7r“7r°/7r°7r°7r° provide a means to extract the 
d, u quark mass difference nid — rriu [10]. 

This paper reports a study of Jjip —>■ (^7r+7r“7r° and 
J /ip —>■ (p'K^'K^'K^ with (p —>■ K^K~ based on a sample of 
(1.311 ± 0.011) X 10® [11, 12] J/r/j events accumulated with 
the BESIII detector in 2009 and 2012. 

II. DETECTOR AND MONTE CARLO SIMULATION 

The BESIII detector [13] is a magnetic spectrometer lo¬ 
cated at the Beijing Electron-Positron Collider (BEPCII), 
which is a double-ring e+e“ collider with a design luminos¬ 
ity of 10^® cm“®s“^ at a center of mass (c.m.) energy of 


* The PDG [1] gives an average value, V(ri' —>■ 37r°)/r(?7' —>■ = 

0.0078 ± 0.0010, of three measurements [7-9]. B{r}' —>■ Stt^) is calcu¬ 
lated using B{ti' —>■ TT^TT^T/) = 0.222 ± 0.008 [1], assuming the uncer¬ 
tainties are independent. 


3.773 GeV. The cylindrical core of the BESIII detector con¬ 
sists of a helium-based main drift chamber (MDC), a plastic 
scintillator time-of-flight system (TOE), and a CsI(Tl) electro¬ 
magnetic calorimeter (EMC). All are enclosed in a supercon¬ 
ducting solenoidal magnet providing a 1.0 T (0.9 T in 2012) 
magnetic field. The solenoid is supported by an octagonal 
flux-return yoke with resistive plate counter muon identifier 
modules interleaved with steel. The acceptance for charged 
tracks and photons is 93% of Ttt solid angle. The charged- 
particle momentum resolution is 0.5% at 1 GeV/c, and the 
specific energy loss (dE/dx) resolution is better than 6%. The 
photon energy is measured in the EMC with a resolution of 
2.5% (5%) at 1 GeV in the barrel (endcaps). The time res¬ 
olution of the TOE is 80 ps (110 ps) in the barrel (endcaps). 
The BESIII offline software system framework, based on the 
Gaudi package [14], provides standard interfaces and utilities 
for event simulation, data processing and physics analysis. 

Monte Carlo (MC) simulation, based on the GEANT4 [15] 
package, is used to simulate the detector response, study the 
background and determine efficiencies. Eor this analysis, we 
use a phase space MC sample to describe the three body decay 
J /ip ^ 077®/0 (980), while the angular distributions are con¬ 
sidered in the decays J/ip —>■ (()/i(1285) ^ (pn^ and 
J/ip —>■ (pT]' . In the MC samples, the width of the /0(98O) is 
fixed to be 15.3 MeV/c®, which is obtained from a fit to data as 
described below. An inclusive MC sample of 1.2 billion J/ip 
decays is used to study the background. Eor this MC sam¬ 
ple, the generator BesEvtGen [16, 17] is used to generate 
the known J/ip decays according to their measured branching 
fractions [1] while LUNDCHARM [18] is used to generate the 
remaining unknown decays. 


III. EVENT SELECTION 

Charged tracks are reconstructed from hits in the MDC 
and selected by requiring that |cos0| < 0.93, where 0 is 
the polar angle measured in the MDC, and that the point 
of closest approach to the e“''e“ interaction point is within 
±10 cm in the beam direction and within 1 cm in the plane 
perpendicular to the beam direction. TOE and dE/dx infor¬ 
mation are combined to calculate the particle identification 
(PID) probabilities for the pion, kaon and proton hypothe¬ 
ses. Eor each photon, the energy deposited in the EMC must 
be at least 25 MeV (50 MeV) in the region of | cos0| < 0.8 
(0.86 < I cos0| < 0.92). To exclude showers that originate 
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from charged tracks, the angle between a photon candidate 
and the closest charged track must be larger than 10°. The 
timing information from the EMC is used to suppress elec¬ 
tronics noise and unrelated energy deposits. 

To be accepted as a J/'0 K'^ decay, a 

candidate event is required to have four charged tracks with 
zero net charge and at least two photons. The two oppositely 
charged tracks with an invariant mass closest to the nominal 
mass of the </> are assigned as being kaons, while the remain¬ 
ing tracks are assigned as being pions. To avoid misidentifi- 
cation, kaon tracks are required to have a PID probability of 
being a kaon that is larger than that of being a pion. A 5- 
constraint kinematic fit is applied to the candidate events un¬ 
der the hypothesis J/'0 A'+iT“7r+7r“77. This includes a 

constraint that the total four-momenta of the selected particles 
must be equal to the initial four-momentum of the colliding 
beams (4-constraint) and that the invariant mass of the two 
photons must be the nominal mass of the 7 r° ( 1 -constraint). 
If more than 2 photon candidates are found in the event, the 
combination with the minimum x^(5C') from the kinematic 
fit is retained. Only events with a x^(5C') less than 100 
are accepted. Events with a K^tt^ invariant mass satisfy¬ 
ing |M(iT^ 7 r^) — M{K*^)\ < 0.050 GeV/c^ are rejected 
in order to suppress the background containing or 
intermediate states. 


To be accepted as a J/'0 Ar+Ar“ 7 r° 7 r° 7 r° decay, a can¬ 

didate event is required to have two oppositely charged tracks 
and at least six photons. Eor both tracks, the PID probability 
of being a kaon must be larger than that of being a pion. The 
six photons are selected and paired by minimizing the quan¬ 


(M (7172) — 


(M(7374) —M„o)^ 


(M (7576)—M„o 


tity 

''-k" ~-k" "ttU 

where is the mass of 7 i 7 j, and M^o and tr^o are 

the nominal mass and reconstruction resolution of the 7 r° re¬ 
spectively. A 7-constraint kinematic fit is performed to the 
J/ij} —?> hypothesis, where the constraints include 

the four-momentum constraint to the four-momentum of the 
colliding beams and three constraints of photon pairs to have 
invariant masses equal to the 7 r°. Events with a x^(7C') less 
than 90 are accepted. 


Eigures 1 (a) and (b) show M{3tt) versus M{K^K~) for 
the two final states respectively. Clear signals from (prj and 
4)7]' with 77 ' — 37 r° are noticeable. In Eig. 1 (a), horizontal 
bands are noticeable from uj and (j) decaying into 7 r+ 7 r“ 7 r° in 
the background channel J /'ip ^ uj / (pK'^ K~. 

To search for the decay J /ip —>■ (/)7r*’/o(980), we focus on 
the region 0.99 < M{K+K-) < 1.06 GeV/c^ and 0.850 < 
< 1.150 GeV/c^. TheM(A'+Ar“) spectra are shown 
in Eig. 2. Clear (p signals are visible. The M( 7 r+ 7 r“) and 
M ( 7 r° 7 r°) spectra for the (p signal region, which is defined by 
requiring 1.015 < M{K~^K~) < 1.025 GeV/c^, are pre¬ 
sented in Eig. 3 (a) and (b) respectively. A clear /o(980) 
peak exists for the tt+tt” mode. The M(/o(980)[7r7r]7r°) 




FIG. 1. Scatter plots of (a) M( 7 r+ 7 r 7 r°) versus M{K^ K ) and 
(b) M{'K^'K°Tr^) versus K~). 


spectra for the /o(980) signal region, defined as 0.960 < 
M{'K'k) < 1.020 GeV/c^, are presented in Eig. 4. There is 
evidence of a resonance around 1.28 GeV/c^ for the decay 
/o(980) —>■ 7 r+ 7 r“, which will be identified as the /i(1285) 
To ensure that the observed /o and /i signals do not orig¬ 
inate from background processes, the same selection criteria 
as described above are applied to an MC sample of 1.2 billion 
inclusive J/ip decays which does not contain the signal de¬ 
cay. As expected, neither an /i nor an /o is observed from the 
inclusive MC sample. The non-c^ background is studied using 
data from the (p sideband regions (0.990 < M{K^K~) < 
1.000 GeV/c^ and 1.040 < M{K+K-) < 1.050 GeV/c^), 
which are given by the hatched histograms in Eig. 3 and Eig. 4 
and in which no /q or /i signals are observed. 


IV. SIGNAL EXTRACTION OF J/ip cpTv°fo (980) 

Eigures 3 (a) and (b) show the tt+tt” and mass spectra 
for events with M{K~^K~) in the (p signal region (the black 


^ For simplicity, /o(980) and /i(1285) will be written as /o and fi respec- 
lively throughout this paper. 
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dots) and sideband regions (the hatched histogram scaled by 
a normalization factor, C). Events in the (p sideband regions 
are normalized in the following way. A fit is performed to 
the K^K~ mass spectrum, where the (p signal is described 
by a Breit-Wigner function convoluted with a Gaussian reso¬ 
lution function and the background is described by a second- 
order polynomial. The mass and width of the p resonance are 
fixed to their world average values [1] and the mass resolu¬ 
tion is allowed to float. The normalization factor C is defined 
as Asig/Asbd, where Asig (Asbd) is the area of the background 
function from the fits in the signal (sideband) region. The re¬ 
sults of the fits are shown in Fig. 2 (a) and (b). 

To extract the signal yield of Jftp pir'^fo, a simulta¬ 

neous unbinned maximum likelihood fit is performed to the 
7r+7r“ and 7r°7r° mass spectra. The lineshape of the /o signal 
is different from that of the Flatte-form resonance observed in 
the decays J/tp pTr'^TT~ and J/p pK'^K~ [19]. A 
Breit-Wigner function convoluted with a Gaussian mass reso¬ 
lution function is used to describe the /o signal. The mass res¬ 
olutions of the fo in the and M(7r°7r°) spectra are 

determined from MC simulations. The non-i/) background is 
parameterized with a straight line, which is determined from 
a fit to the data in the p sideband regions. The size of this 
polynomial is fixed according to the normalized number of 
background events under the p peak, Abkg = GWsbd, where 
Nshd is the number of events falling in the p sideband regions 
and C is the normalization factor obtained above. Another 
straight line is used to account for the remaining background 
from J/p ^ pTr^TiTT without /o decaying into tttt. 

The mass and width of the /o are constrained to be the same 
for both the and the Ar+Ar“7r°7r°7r° final 

states. The fit yields the values M(/o) = 989.4 ± 1.3 MeV/c^ 
and r(/o) = 15.3 ± 4.7 MeV/c^, with the number of events 
N = 354.7 ±63.3 for the tt+tt” mode and 69.8 ±21.1 for the 
7r°7r° mode. The statistical significance is determined by the 
changes of the log likelihood value and the number of degrees 
of freedom in the fit with and without the signal [20]. The 
significance of the /o signal is 9.4(T in the 
final state and 3.2 (t in the Ar+iT“7r°7r°7r° final state. The 
measured mass and width obtained from the invariant dipion 
mass spectrum are consistent with those from the study of the 
decay J/p —>■ 7?7(1405) ^ 77r°/o(980) [4]. It is worth not¬ 
ing that the measured width of the /o observed in the dipion 
mass spectrum is much smaller than the world average value 
of 40-100MeV [1]. 


V. SIGNAL EXTRACTION OF J/p pfi (1285) WITH 
/i(1285) ->■ 7r°/o(980) 

Figures 4 (a) and (b) show the 7r+7r“7r° and mass 

spectra in the p and /q signal region (the black dots) and side¬ 
band regions (the hatched histogram). The /o sideband re- 




FIG. 2. Fits to the M(A'+A") mass spectra for the mode (a) 
/o(980) —>■ 7r''"7r“ and (b) /o(980) —^ The solid curve is the 

full fit; the long-dashed curve is the p signal while the short-dashed 
curve is the background. 

gions are defined as 0.850 < M{TnT) < 0.910 GeV/c^ and 
1.070 < Mpm) < 1.130 GeV/c^. In Fig. 4, events in the 2- 
dimensional sideband regions are weighted as follows. Events 
that fall in only the p or /o(980) sideband regions are given 
a weight 0.5 to take into account the non-(() or non-/o(980) 
background while those that fall in both the p and the /o(980) 
sideband regions are given a weight —0.25 to compensate 
for the double counting of the non-t/) and non-/o(980) back¬ 
ground. There is evidence of a resonance around 1.28 GeV/c^ 
that is not noticeable in the 2-dimensional sideband regions. 
By studying an MC sample of the decay J/p pfi 
anything, we find that the decay /i Tr^Tr^p/Tr^ao ^ with 
■7 77 contributes as a peaking background for the decay 

/i 7r°7r°7r°. The yield of this peaking background is calcu¬ 
lated to be 3.1 ±0.6 using the relevant branching fractions"^ [1] 
and the efficiency determined from an MC simulation. A si¬ 
multaneous unbinned maximum likelihood fit is performed to 


^ For simplicity, ag (980) and a[| (980) are written as oq and respectively 
throughout this paper. 

'* We assume that B{fi —^ Tr^Tr^rj) = ^B(fi —>■ Trwr/), B{fi —)■ 
Tr'^a^) = ^B{fi —> Trao), and B(a§ —>■ ir^rj) = 100%. 
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FIG. 3. The spectra (a) M( 7 r+ 7 r“) and (b) M( 7 r° 7 r°) (three entries 
per event) with K^K~ in the cj) signal region (the black dots) and 
in the 0 sideband regions (the hatched histogram). The solid curve 
is the full fit; the long-dashed curve is the /o(980) signal; the dotted 
line is the non-<j) background and the short-dashed line is the total 
background. 

the M(7r+7r“7r°) and M(7r°7r°7r°) distributions. The /i sig¬ 
nal is described by a Breit-Wigner function convoluted with a 
Gaussian mass resolution function. The shape of the peaking 
background fi r ]/is determined from an exclu¬ 

sive MC sample and its size is fixed to be 3.1. A second order 
polynomial function is used to describe the remaining back¬ 
ground. The mass resolutions of the /i in M(7r+7r“7r°) and 
M(7r°7r°7r°) are determined from MC simulations. 

The fit to M(7r+7r“7r°) and M(7r°7r°7r*i) distributions 
yields the values M[fi) = 1287.4 ± 3.0 MeV/c^ and 
r(/i) = 18.3 ± 6.3 MeV/c^, with the number of events 
N = 78.2 ± 19.3 for the iT+/C“7r+7r“7r° final state and 
(V = 8.7 ± 6.8 (< 18.2 at the 90% Confidence Level (C. 
L.)) for the iC+iT“7r°7r°7r° final state. The mass and width 
are consistent with those of the axial-vector meson fi [1] 
The statistical significance of the fi signal is 5.2cr for the 


^ Here we assume that the contribution of the pseudoscalar r)(1295) is small 
as no significant r7(1295) signals were found in the 7r+7r~r; mass spectrum 
from a study of J/tp —> <)>7r+7r“»7 [21]. 


FIG. 4. The spectra of (a) and (b) M( 7 r° 7 r° 7 r°) in the 

4> and /o (980) signal region (the black dots with error bars) and in 
the sideband regions (the hatched histogram). The solid curve is the 
result of the fit, the long-dashed curve is the /i(1285) signal, and 
the short-dashed curve is the background. In (b), the dotted curve 
represents the peaking background from the decay /i(1285) —^ 
Tr'^TT^Ty/TT^ao with 77 —>■ 77 . 

K'^K~Tr~^'K~Tr'^ final state and 1.8 (t for the K~ 
final state. From the fit results, summarized in Table I, it is 
clear that the production of a single fi resonance cannot ac¬ 
count for all of the /ott® events above the background. 

VI. SIGNAL EXTRACTION OF J/tp rprj' 

For the decay J/ij) cf>r]' — K'^K~TT~^Tr~Tr'^, the de¬ 
cays J/ij) —>■ 4)7]' —?> and J/ijj —>■ 

K~'K~ produce peaking background. 
To reduce the former peaking background which is dominant, 
events with 0.920 < < 0.970 GeV/c^ are re¬ 

jected. 

As the amount of background for the decay J/tjj ^ (f’V' 
K'^K~7T^Tr^Tr^ is relatively small, the (f) signal and side¬ 
band regions are expanded to be 1.010 < M{K+K-) < 
1.030 GeV/c^ and 1.040 < M{K+K-) < 1.060 GeV/c^, re¬ 
spectively. A peaking background for this decay comes from 
the decay J/'0 (pr]' —>• To reduce this 
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TABLE I. Summary of the observed number of events the 

errors are statistical only.). 


obs 


Decay mode 


N' 


fo 

07r°/o, fo 


J/i’- 

J/i’- 


j/tA ^ 0/1,/i 


n°fo,fo - 
^°fo, fo 


354.7 ± 63.3 
69.8 ± 21.1 
78.2 ± 19.3 

8.7 ± 6.8 

< 18.2 (90% C.L.) 


+ = cr 

TT TT TT 


j/if 4 >v'’ v' 

J/ip cl)ri', Tj' —>■ 7r°7r°7r° 


183.3 ± 21.0 
77.6 ±9.6 



background, events with any photon pair mass in the range 
0.510 < M( 77 ) < 0.580 GeV/c^ are rejected. 

Figures 5 (a) and (b) show the final 7 r+ 7 r“ 7 r° and 
mass spectra for the (j) signal (the black dots) and sideband 
(the hatched histogram) regions. By analyzing data in the (f 
sideband regions and the inclusive MC sample, we find that 
the contribution from the decay J/ip K~^K~t]' is negligi¬ 
ble. 

An unbinned likelihood fit is performed to obtain the sig¬ 
nal yields. The rj' signal shape is determined by sampling a 
histogram from an MC simulation convoluted with a Gaus¬ 
sian function to compensate for the resolution difference be¬ 
tween the data and the MC sample. The shape of the peak¬ 
ing background is determined from exclusive MC samples, 
where the relative size of the background shape is determined 
using the relevant branching fractions in the PDG [1]. The 
non-peaking background is described by a first-order (zeroth- 
order) polynomial for the 77 ' 7 r+ 7 r“ 7 r° decay. 

The number of events are determined to be TV = 183.3 ± 21.0 
for the iC+TV“ 7 r+ 7 r“ 7 r*’ final state and 77.6 ± 9.6 for the 
K'^K~TT^TT^Tr^ final state. 



FIG. 5. The spectra (a) M{'k^'k 7 r°) and (b) M( 7 r° 7 r° 7 r°) with 
K+K- in the cj> signal region (the black dots) and sideband regions 
(the hatched histogram). The solid curve is the result of the fit, the 
long-dashed curve is the rj' signal, and the short-dashed line is the 
polynomial background. In (a), the dotted and dot-dashed curves 
represent the peaking background 7 ' —>■ 7 p ^ 7 ( 7 ) 77 ^ 77 “ and 
rj' —>■ 70 ; —^ 777 ^ 71 “TT®, respectively. In (b), the dotted curve repre¬ 
sents the peaking background rj' —>■ with 7 —>■ 77. 


VII. BRANCHING FRACTIONS MEASUREMENT 


B{rj' ^X) = 


T^obs 






( 2 ) 


Table I summarizes the signal yields extracted from the fits 
for each decay. Equations (1) and (2) give the formulae used 
to calculate the branching fractions, where n is the number of 
7 r°s in the final state X. TV°*”* and e are the signal yield from 
the fits and efficiency from the MC simulation for each decay, 
respectively. Byz is the branching fraction of the decay X 
YZ. Nj/^ is the number of J/ijj events. The upper limit 
of —>■ (ffi,/! 7r°/o,/o —f 7 r° 7 r°) is determined 

according to Eq. (3), where TV°pp is the signal yield at the 90% 
C. L. and is the total systematic uncertainty, which is 
described in the next section. Equation (4) is used to calculate 
the ratio between the branching fraction for 7 ' and 

that for 7 ' —> 7 r+ 7 r“ 7 r°. 


BiJ/if (jrX) 


ATobs 


( 1 ) 


B{J/if <j)X) < 


ATOhs 

^^upp 


(3) 


tstt = 6 ( 7 '—>■ 7 r° 7 r° 7 r°)/S( 7 '—TT+TT 7 r°) 
TV°'’'*(7r°7r°7rO) e(7r+7r-7r0) 1 

TV°bs(7r+7r“7rO) e(7rb7r*^7r*^) {BJ^°)'^ 


VIII. ESTIMATION OF THE SYSTEMATIC 
UNCERTAINTIES 


(1) MDC tracking: The tracking efficiency of kaon tracks 
is studied using a high purity sample of J/tp —>■ KsKtt 
events. The tracking efficiency of the low-momentum 
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pion tracks is studied using a sample of J/tp 
'K~^TT~pp while that of the high-momentum pion tracks 
is studied using a high statistics sample of J/ip —>■ pir. 
The MC samples and data agree within 1% for each 
kaon or pion track. 

(2) Photon detection: The photon detection efficiency is 
studied using a sample of J/tp pir events. The sys¬ 
tematic uncertainty for each photon is 1% [22]. 

(3) PID efficiency: To study the PID efficiency for kaon 

tracks, we select a clean sample of J/ip (pp 

K+K- ' 77 . The PID efficiency is the ratio of the num¬ 
ber of events with and without the PID requirement for 
both kaon tracks. MC simulation is found to agree with 
data within 0.5%. 

(4) Kinematic fit: The performance of the kinematic 

fit is studied using a sample J/ip —>■ cpp —> 
iT+iC“7r+7r“7r°/iT+iT“7r°7r°7r°, which has the same 
final states as the signal channel J/ip —)• (pir'^fo with 
(p — K'^K~ and ffi /tt^tt^. The control sam¬ 

ple is selected without using the kinematic constraints. 
We then apply the same kinematic constraints and the 
same requirement on the from the kinematic fit. The 
efficiency is the ratio of the yields with and without the 
kinematic fit. It contributes a systematic uncertainty of 
1.0% for /o —>■ 7r+7r“ and 2.0% for ffi —>• 

(5) Veto neutral K*: In selecting the candidate events 
J/ip —>• (pTT^fo —>■ K^K~Tr~^7r~TT^, the events with 
\M{K^n^)-M{K*°)\ < 0.050 GeV/c^ are vetoed to 
suppress the background containing or K*^ inter¬ 
mediate states. The requirement is investigated using a 
clean sample J /ip <pp ^ A"+A"“7r+7r“7r°. The ef¬ 
ficiency is given by the yield ratio with and without the 
requirement |M(A'^7rT) —M(Ar*°)I < 0.050 GeV/c^. 
The efficiency difference between data and MC is 0.1 %. 

(6) (p signal region: The uncertainty due to the restriction 
on the (p signal region is studied with a high purity sam¬ 
ple of J/ip —?► pp' —s- K'^K~Tr~^Tr~p events as this 
sample is free of the background J/ip K^K~p' 
without the intermediate state p. 

(7) Veto peaking background: The uncertainties due to the 
restrictions used to remove peaking background in the 
mode ry' Stt are studied with a control sample of 
J/ip —>■ ojp —>■ 2(7r+7r“7r°) events. For each sample, 
the efficiency is estimated by comparing the yields with 
and without the corresponding requirement. The differ¬ 
ence in efficiency between the data and MC samples is 
taken as the systematic uncertainty. 

(8) Background shape: To study the effect of the back¬ 
ground shape, the fits are repeated with a different fit 


range or polynomial order. The largest difference in 
signal yield is taken as the systematic uncertainty. 

(9) Mass resolution: The mass resolutions, ctmCj from an 
MC simulation of the modes /o —J- 7r+7r“/7r°7r° and 
/i —>■ 7r°/o have an associated systematic uncertainty. 
The difference in mass resolution, era, between the data 
and the MC simulation is determined using a sample of 
J/ip pp events where p —>■ 7r+7r“7r°/7r°7r°7r°. The 
fit is repeated using different mass resolutions, which 
are defined as assuming ac is the same 

for the two-pion and three-pion mass spectra. The dif¬ 
ference in yield is taken as a systematic uncertainty. 

(10) MC simulation: For the decay J/ip pir'^fo, the 
dominant systematic uncertainty is from the efficiency 
eo determined by a phase space MC simulation. The 
7r°/o invariant mass spectrum is divided into 5 bins, 
each with a bin width of 0.2 GeV/c^. The /q sig¬ 
nal yields, Ni, are determined by fits to the tttt spec¬ 
tra for each bin i using the mass and width of the 
/o obtained above. The corrected efficiency is cm = 

, where Ci is the efficiency in the z-th bin. The 
same procedure is applied to the angular distribution 
of the 7r°/o system in the c.m. frame of the J/ip to 
obtain another corrected efficiency ee- The difference 
{cm — eo)^ + (ee — eo)^ is taken as the systematic 
uncertainty due to the imperfection of the MC simu¬ 
lation. 

(11) /o signal region: For the decay J/ip pfi with 

fi —>■ 7r°/o, the /o signal region is 0.960 < M{'K'k) < 
1.020 GeV/c^. The branching fraction measurements 
are repeated after varying this region to 0.970 < 
Mijr-K) < 1.010 GeV/c^ and 0.950 < Mijn:) < 
1.030 GeV/c^. The differences from the nominal re¬ 
sults are taken as the systematic uncertainties due to the 
signal region of the /q. For the decay /i —>■ the 

number of the peaking background fi —^ 7r°7r°p[77] is 
determined to be 3.1 ± 0.6. Varying the number of the 
peaking background within ±0.6 in the fit, the largest 
difference of the signal yield gives a systematic uncer¬ 
tainty. The systematic uncertainty values related to the 
fi are shown in brackets in Table II. 

(12) About B{J/ip PhJi 7r°7r°): For 

the decay J/ip -)> pfiji 7r°/o with ffi 7r°7r°, 
the signal yield at the 90% C. L., A^°pp in Eq. (3), is 
the largest one among the cases with varying the fit 
ranges, the order of the polynomial describing the back¬ 
ground, the number of the peaking background, and 
the signal region of the /o resonance. The total sys¬ 
tematic uncertainty, in Eq. (3), is the quadratic 
sum of the rest systematic uncertainties in the third col¬ 
umn of Table II (the values in the brackets). We ob- 
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tain 7V°pp = 29.0 and = 6.9% with the efficiency 
(7.21 ± 0.08)%, determined from an MC simulation. 

7r°/o,/o 7r°7r°) is calculated 

to be less than 6.98 x 10“^ at the 90% C. L. according 
to Eq. (3). 

(13) Uncertainty of B{J / %!) —>• fv')'- For the decay rj' —>■ Stt, 

the dominant systematic uncertainty arises from the un¬ 
certainty of B{J/'ip — >■ frj') = (4.0 ± 0.7) x 10“'^ [1]. 
A variation in S(J/'0 frj') will change the size 
of peaking background and thus the signal yield. In 
Eq. (2), it is reasonable to consider a change in the quan¬ 
tity with any variation in B{J/'tjj ft]'). 

The fit to the data is repeated after varying the num¬ 
ber of peaking background to correspond with Icr vari¬ 
ations in B{J/'ip — >■ frj') [1]. The largest difference 
of from the nominal result is taken as the 

systematic uncertainty. 

(14) Systematic uncertainties for In the measure¬ 

ment of the ratio of B{r]' —>■ tt^tt^tt^) over 
B{r]' —7r+7r“7r°), the systematic uncertainties due 
to the reconstruction and identification of kaon tracks 
and photon detection cancel as the efficiency ratio 
e(7r°7r°7r°)/e(7r+7r“7r°) appears in Eq. (4). The effect 
of the uncertainty in the number of peaking background 
due to the uncertainty of B{J/fi frj') is also consid¬ 

ered. 

All systematic uncertainties including those on the number of 
J/ij} events [12] and other relevant branching fractions from 
the PDG [1] are summarized in Table II, where the total sys¬ 
tematic uncertainty is the quadratic sum of the individual con¬ 
tributions, assuming they are independent. Efficiency and 
branching fraction measurements are summarized in Table III. 

IX. SUMMARY 

In summary, we have studied the decay J/tjj ^ fSir 
K'^K~2>'k. The isospin violating decay J/ft 
is observed for the first time. In the tt^/q mass spec¬ 
trum, there is an evidence of the axial-vector meson /i, 
but not all 7r°/o pairs come from the decay of an /i. Us¬ 
ing BiJ/ilt ffi) = (2.6 ± 0.5) X 10-4 and B{fi 

TTUo —>■ TTTrry) = (36 ± 7)% from the PDG [1], the ratio 
B{fi 7 r°/o —?► 7r°7r+7r-)/B(/i —Tr^ag — n^TT^r]) is 
determined to be (3.6 ± 1.4)% assuming isospin symmetry 
in the decay fi —>■ aoir. This value is only about 1/5 of 
^(77(1405) 7r°/o —>■ 7r°7r+7r-)/;B(r7(1405) tt^Oq 

TT^TT^ri) = (17.9±4.2)% [4]. On the other hand, the measured 
mass and width of the /o obtained from the invariant dipion 
mass spectrum are consistent with those in the decay J/ip ^ 


777 ( 1405 ) —>■ 77r°/o [4]. The measured /g width is much 
narrower than the world average value of 40 — 100 MeV [1]. 
It seems that there is a contradiction in the isospin-violating 
decays /i/77(1405) —>■ 7r°/g. However, a recent theoretical 
work [23], based on the triangle singularity mechanism as 
proposed in Ref. [5, 6], analyzes the decay fi tt^fo —t 
7r°7r+7r- and predicts that the width of the peaking structure 
in the /g region is about 10 MeV. It also derives B{fi —)■ 
7r°/g —>• 7r°7r+7r-)/S(/i —>■ Tr^Ug —>■ n^n^r]) ~ 1%, which 
is close to our measurement. This analysis supports the argu¬ 
ment that the nature of the resonances Og and /g as dynami¬ 
cally generated makes the amount of isospin breaking strongly 
dependent on the physical process [23]. In addition, we 
have measured the branching fractions B{ri' tt+tt-tt^) = 
(4.28 ± 0.49(stat.) ± 0.22(syst.) ± 1.09) x lO-^ and B{i 
ttOttOttO) = (4.79 ± 0.59(stat.) ± 0.33(syst.) ± 1.09) x IQ-^, 
where the last uncertainty is due to B{J/ ip —?> prj'). The ratio 
between them = 1.12 ± 0.19(stat.) ± 0.06(syst.) is also 
measured for the first time. These results are consistent with 
those measured in the decay J/ip ^ 777 ' [4]. 
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TABLE II. Summary of systematic uncertainties (%). For /o ^ tttt, the values in the brackets are for the decay /i ^ T^^fo- For t]' —>■ Stt, the 
systematic uncertainty from the uncertainty of B{J/ip 4>'n') is not included in the total quadratic sum. The last column lists the systematic 
uncertainties for the ratio between B{rj' —>■ 7 r° 7 r'^ 7 r°) and B{rj' —>■ denoted by 


Sources 

fo tt+tt” 

fo 7 r®7r® rj' 

— ( 

—¥ 7T TV TT 

T]' StT^ 

ran 

MDC tracking 

4.0 

2.0 

4.0 

2.0 

2.0 

Photon detection 

2.0 

6.0 

2.0 

6.0 

4.0 

PID efficiency 

0.5 

0.5 

0.5 

0.5 

- 

Kinematic fit 

1.0 

2.0 

1.0 

2.0 

1.5 

Veto neutral K* 

0.1 

- 

- 

- 

- 

(p signal region 

1.1 

1.1 

1.1 

0.5 

0.5 

Veto peaking bkg. 

- 

- 

0.3 

0.9 

0.9 

Bkg. shape 

5.4(15.5) 

4.4(15.6) 

1.3 

0.3 

1.4 

Mass resolution 

0.3 (0.4) 

1.0 (0.1) 

- 

- 

- 

MC simulation 

11.4 (-) 

11.4 (-) 

- 

- 

- 

fo signal region 

-(2.4) 

-(68.2) 

- 

- 

- 

B{j/p, ^ tpn') 

- 

- 

25.6 

22.8 

- 

Peaking bkg. 

- 

-(6.9) 

- 

- 

2.2 

Number of J/ip 

0.8 

0.8 

0.8 

0.8 

- 

Other B.F. 

1.0 

1.0 

1.0 

1.0 

0.1 

Total 

13.6(16.5) 

14.5 (70.6) 

5.1 

6.9 

5.5 


TABLE III. Summary of the efficiencies and the branching fractions. Eor the branching fractions, the first error indicates the statistical error 
and the second the systematic error. For B{ri' ^ Stt), the third error is due to the uncertainty of B{J/ip —)■ tprj') [ 1 ]. The last line gives the 
measured value of r3,r, defined as B{t]' —>■ 'K°Tv^Tr^)/B{T]' —^ 7r+7r“7r°). 


Decay mode Efficiency (%) Branching fractions 

J/ip (pn'^fojo tt+tt” 12.44 ± 0.10 (4.50 ± 0.80 ± 0.61) x 10"'^ 

J/ip (pn°fo, fo n°Tv° 6.76 ± 0.08 (1.67 ± 0.50 ± 0.24) x 10"® 

J/ip (pfi, /i TT°fo 7r°7r+7r- 13.19 ± 0.11 (9.36 ± 2.31 ± 1.54) x lO"’’ 

J/ip (pfi, fi TT°fo TT°TT°n° 6.76 ± 0.08 (2.08 ± 1.63 ± 1.47) x 10"’^ 

< 6.98 X 10"'^ (90% C. L.) 


rj' 7V~^7V 7V^ 

rj' 7r®7r®7r® 

16.92 ± 0.12 (4.28 ± 0.49 ± 0.22 ± 1.09) x 10"® 
6.55 ± 0.08 (4.79 ± 0.59 ± 0.33 ± 1.09) x 10"® 

ran 

1.12 ±0.19 ±0.06 
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